function c = p1(t)
    c = [1 -1 0 t^3];